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Abstract. In this work the phenomenology of models possessing a non-minimal coupling 
between matter and geometry is discussed, with a particular focus on the possibility of describing 
the flattening of the galactic rotation curves as a dynamically generated effect derived from this 
modification to General Relativity. Two possibilities are discussed: firstly, that the observed 
discrepancy between the measured rotation velocity and the classical prediction is due to a 
deviation from geodesic motion, due to a non-(covariant) conservation of the energy-momentum 
tensor; secondly, that even if the principle of energy conservation holds, the dynamical effects 
arising due to the non-trivial terms in the Einstein equations of motion can give rise to an 
extra density contribution that may be interpreted as dark matter. The mechanism of the 
latter alternative is detailed, and a numerical session ascertaining the order of magnitude of the 
relevant parameters is undertaken, with possible cosmological implications discussed. 



1. Introduction 

Many modifications of the tlieories of gravity are strongly motivated by the issue of an alternative 
solution to the problem of the flattening of the galaxy rotation curves — one which does not 
require dark matter in the strictest sense, i.e. weakly or non-interacting matter fields [Tl[2l[3lll]- 
With this goal in mind, standard f{R) models [5l[6l[7] have been studied endowed with power- 
law curvature terms /i(-R) oc i?" in the action, instead of the linear curvature depicted in 
the General Relativity (GR) action. For n > 1, an addition to the Newtonian potential 
A$(r) = —GM{r /tcY /2r is found, with /? a function of n and Vc being a parameter characteristic 
of each galaxy. A fit of several rotation curves indicates that n = 3.5 (corresponding to 
/3 = 0.817) yields the best agreement with observations [8j. 

In another approach, it has been shown that an asymptotically flat velocity dispersion curve 
can be obtained if the curvature term in the Einstein-Hilbert action includes a logarithmic factor, 
fi{R) = /o-R(l + log i?) j9j; given the non-relativistic condition v <^1, this may be approached 
by a power-law with an exponent very close to unity, fi{R) ~ foR^^^ . 

Since this description does not take into consideration any galaxy-dependent parameters, an 
universal asymptotic velocity v is obtained; this is in stark contrast with the Tully-Fisher and 
Faber-Jackson relations, which posit a dependence of the visible mass of a galaxy M oc f™, for 
spiral (m = 4) and elliptical galaxies (m = 6), respectively. 



^ Also at Departamento de Fi'sica, Instituto Superior Tecnico. 



This work reviews the results obtained in Ref. [10], to which the reader is directed for a 
more technical discussion and the full set of computations. It comprises a small introduction to 
the non-minimal coupling model under scrutiny, proceeds with the implied non-geodesic motion 
already hinted above, and then proposes the already mentioned mechanism for dynamically 
mimicking dark matter; a numerical session and discussion follows. 



2. The model 

The non-minimal coupling between matter and geometry is embodied in the following action 
functional 



S 



fl{R) + [1 + Xf2iR)]Cr. 
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(1) 



where fi{R) (with i = 1,2) are arbitrary functions of the scalar curvature i?, Cm is the 
Lagrangian density of matter and g is the metric determinant. The non-minimal coupling 
/2 is gauged via the coupling constant A (with dimensions [A] = [/2]~^). The standard Einstein- 
Hilbert action is recovered by taking /2 = and /i = 2k,[R — 2A), where k = c^/lGvrG and A is 
the cosmological constant (from now on, one works in a unit system where c = 1). 

The modified Einstein equations of motion arise from variation of the action with respect to 
the metric g^^, 
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defining \3^u = V^V,y for convenience, and Fi{R) = fi{R). The matter energy- momentum 
tensor is defined as 
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so that the trace of Eq. ([2]) reads 



(Fi + 2AF2£„) R-2h= - 3n (Fi + 2AF2£„) + (1 + A/2) T. 



The Bianchi identities, V^G 
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imply the covariant "non-conservation law" 
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interpreted as a signal of energy-momentum exchange between matter and geometry due to the 
non-trivial terms [12J . 

For simplicity, and since one aims studying the phenomenology of the non-minimal coupling 
f2{R), one sets fi{R) = 2kR (the GR form discarding the cosmological constant A); Eq. ([2]) 
then reduces to 
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and, equivalently, Eq. (jH) becomes 
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3. Flattening of the galaxy rotation curves due to non-geodesical motion 

Prom a fundamental standpoint, the most striking cliaracteristic of ttie model here considered 
is the non-conservation of the energy- momentum tensor of matter; assuming that matter is 
described by a perfect fluid distribution with the usual energy-momentum tensor 



{p+p)UfJJy+Pgf^y, 



(8) 



where p is the energy density, p{p) is the pressure (given by a suitable equation of state) and 
Ufj_ is the four- velocity (satisfying the normalization condition U^W^ = —1), Eq. dSJ leads to an 
extra force imparted upon a test particle. 
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Qfiu — U^Uu, so that hf^yU^ = and the extra force is 



with the projection operator h^i, 
orthogonal to the four-velocity. 

In GR, the specific form of the Lagrangian density does not appear explicitly in the equations 
of motion, but only via the derived energy- momentum tensor Eq. For this reason, and 

particularly if one is not interested in describing the thermodynamics of the perfect fluid, this 
functional is usually ignored. This, of course, changes in the present situation, since Cm appears 
explicitly in Eqs. ^ and Q. 

In GR, the Lagrangian density Cm of a perfect fluid several equivalent on-shell expressions 
are admitted: for instance, Cmo = —p-, Cmi = P and Cm2 = —na, where n is the particle number 
density, and a is the physical free energy, defined as a = p/n — Ts, with T the fluid temperature 
and s the entropy per particle. However, in the presence of a non-minimal gravitational coupling 
with matter, this degeneracy is lifted and the original, bare Lagrangian density Cm = —p should 
be considered il3][Ti|. 

One may further simplify computations by noticing that as matter in galaxies is non- 
relativistic (asides from the inner core region), and it may be described approximately by a 
pressureless dust distribution, p <C pc^- 

With the above consideration in mind, one now assumes a static, spherically or cylindrically 
symmetric Ansatz for the problem and that the velocity has only a tangential component, so 
that h^^ = g^^ and so the extra-force Eq. ([9]) reads 
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where the prime denotes differentiation with respect to the radial coordinate r. 

Endowed with the expression above, one may now t le the issue of describing the flattening of 
the galaxy rotation curves; instead of an asymptotic velocity dispersion v{r) — )■ Uqo, one assumes 
that the reported flattening allows for a shallow decrease, i.e., a profile of the form 



Vr^ — 
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where r* is an unknown normalization constant and < /3 < 1 is the outer slope of the rotation 
curve. Identifying the extra- force with the observed deviation leads to the expression 
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which one may use to ascertain the form of the non-minimal coupling f2{R)', in order to do 
so, one requires some bearing on the dependence of the scalar curvature on the radial distance 
R = R{r); to do so, one recalls that it is implicitly assumed that the effect of the non-trivial 



terms appearing in Eq. ([6|) is perturbative — so that the geodesies derived from the metric are 
similar to those found in GR — and the non-minimal coupling manifests itself simply through 
an extra-force deviating test particles from geodesic motion. 

As will be seen below, this approach is self-consistent; for now, one assumes that the scalar 
curvature is given by its classical expression R = —SirGT, with T = —p the trace of the energy- 
momentum tensor T^^. Thus, the issue of obtaining a dependence R{r) simplifies to ascribing 
a particular profile for the density distribution p(r) of matter. For this, one takes a power-law 
profile of the form 



p(r) po{a/ry 



(13) 



with m 7^ and po, a characteristic scales. 

With the differential equation Eq. (|12p in a closed form, one may solve for f2if') and, reversing 
the dependences above, find the aimed coupling f2{R)- After some algebraic steps (detailed in 
Ref. [lO]), one finds that 
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defining the characteristic scale R^, = 8'irGpo{a/r^:)"^. 

With the above solution, it is trivial to show that the conditions for a negligible effect of the 
non-minimal coupling f2{R) on the Einstein field Eq. ([6]), 



|A/2(i?)| « 1 
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(using that T = Cm = —p) are satisfied for a wide domain, 1 kpc ^ r < kpc, where the 
fiattening of the rotation curves is measured for a wide range of values of (mainly due to the 
smallness of the gauging constant A = v'^ « 10~^^): one could set, e.g., r^, = 1 kpc. 

As a result, the assumption of a perturbative regime for the equations of motion Eq. ([6]) is 
indeed self-consistent — although the deviation from geodesical motion is relevant, dominating 
the Newtonian acceleration and leading to the flattening of the galaxy rotation curve. 

Before continuing, one remarks that a similar exercise would ensue if one was to assume that 
the velocity dispersion curve follows a logarithmically decaying proflle. 



V [r) = Woo— log — 
r V r* 
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instead of the simpler power-law assumed above. Indeed, and in the context of standard galactic 
dark matter descriptions, this radial dependence leads to a dark matter density distribution 
given by the so-called Navarro-Frenk- White (NEW) density profile, which is widely used in 
the literature |15j (see Appendix A of Ref. [lOj for a discussion of different density profiles). 
Assuming this decay law for f^(r) leads to the non-minimal coupling 



l/m 



1 + — log — 
m \ R 



(17) 



and, as before, the perturbative conditions Eq. ()15p are trivially satisfied within the relevant 
range. 

To close this section, one first remarks that the obtained expression Eq. ()14p for the non- 
minimal coupling f2{,R) in the fiattened rotation curve case /3 = is strikingly similar to that 



of Ref. [9], where the non-trivial curvature term is given by fi{R) = R[l + v'^ log{R / R^^)] w 

2 

R{R/ R^Y°° . Following that proposal, one may suggestively rewrite the action Eq. ([T]) as 



defining a = v^jm. 

However, this similarity with Ref. [U] leads to the same caveat, namely an undesired 
universality of the asymptotic value of the predicted galaxy rotation curves: assuming similar 
outer slopes m for the visible matter density profile p of galaxies, one concludes that v^o and /3 
should be fixed parameters of the model, unadjustable for different galaxies. This is not observed 
for the available rotation curves, in which the asymptotic value Voo may vary somewhat (although 
of the order of magnitude v^o ~ lOO/cm/s). One may assume that for different density profiles, 
spatial symmetry and morphology considerations would contribute to alleviate this issue, but 
the obtained universality remains a noticeable disadvantage of this approach — thus prompting 
the search for alternatives within the ample phenomenological arsenal provided by the model 
under scrutiny. 

4. Pressureless dust with non-minimal gravitational coupling 

Motivated by the considerations put forward in the last paragraph, one now assumes the opposite 
stance: instead of considering a perturbative effect of the non-minimal coupling, but a large 
deviation from geodesic motion, one posits that the energy-momentum tensor of matter is 
approximately conserved, but the geodesies thus followed by test particles differ significantly 
from those of GR. As in the previous section, this approach also proves to be self- consistent. 

Given the above, the problem becomes similar to those faced by standard /(-R) theories (that 
is, with f2{R) = 1): solve the Einstein field equations so that the obtained metric gives rise to 
the observed galaxy rotation curves. Since one ultimately aims to describe dark matter as a 
dynamical effect stemming from the action functional Eq. ([1]), it is productive to try to express 
the effect of the non- minimal coupling in terms of an equivalent, mimicked dark matter density. 
This also allows for a direct comparison with existing profiles for the distribution of dark matter, 
arising both from simulations as well as observational evidence. 

For generality, a power law coupling of the form 



is assumed, with n a yet unspecified exponent and Rq a characteristic scale (which amounts to 
a rewriting of the gauging constant A). 

It should be highlighted that Eq. ()19p does not exclude a more elaborate form for the non- 
minimal coupling, which could then be written as a series 



summed over the rational exponents p of the "full" model (integer or not); if this is the case, 
the choice of Eq. ()19p is to be viewed as a valid approximation within the galactic context 
(that is, the dominant term of the series for the corresponding range of R). For this reason, this 
assumption is not in contradiction with the previous work by the authors [16] , which approached 
the possibility of including a linear coupling /2 = R/Ri (assumed relevant in a high density 
environment such as the Sun). 

Again, one models matter within a galaxy as a pressureless perfect fluid, i.e., dust, 
characterized by a density p; the corresponding energy-momentum tensor thus reads 
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Using the Lagrangian density Cr, 
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and Eq. ([7]) becomes 
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As stated at the beginning of this section, the additional terms induced by the non-minimal 
gravitational coupling f2{R) are interpreted as contributing to the dark matter halo density 
profile pdm- This can be better ascertained by ignoring the matter contribution arising from 
GR, so that the non-minimal terms dominate far away the galactic core (where the rotation 
curve flattens). In this exterior region, one may rewrite Eq. ()23p as 
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^.1. Static solution 

Assuming a known visible matter density /9, solving the above equation amounts to obtaining 
the relation R = R{p); this is somewhat analog to the previous exercise, where one assumed 
p{r) and computed R{r). Naturally, the classical identification R = 2Kp fails, since it would 
lead to a differential equation for the visible matter density which might not be satisfied by the 
considered profile p{r). Inspection of Eq. (p^ shows that a solution defined implicitly by the 
expression 
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is self-consistent, since the gradient term vanishes everywhere; accordingly, this is dubbed as a 
"static" solution. This yields 
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introducing the characteristic density po = Rq/2k. 

In order to interpret the additional terms shown in the Einstein field Eqs. ()22|) . one must 
first show that the above "static" solution provides a suitable description of dark matter, i.e. 
the full tensorial aspect of the problem must be taken into account, not merely the algebraic 
identification stemming from the trace Eq. (|24p . Indeed, replacing Eq. ()26p in Eq. ()22p 
allows one to conclude that the dynamic effect due to the non-minimal coupling behaves as a 
perfect fluid characterized by an energy- momentum tensor analog to that of matter, so that this 
mimicked "dark matter" shares the same 4-velocity as visible matter, but has a pressure and 
density given by 
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thus yielding a non-vanishing pressure and an equation of state of the form 
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Pdm = Pdm- (28) 
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Hence, at a sufficiently large distance from the galactic core, where the rotation curve has 
flattened, the non-minimal gravitational coupling with ordinary matter mimics a dark matter 
component which is dragged by ordinary matter, and may have a negative density pdm < 0, if 
1/4 < n < 1; although strange, this is not a pathology of the model, as the curvature in the 
considered outer region will remain positive, i? > 0. As will be seen below, relevant models for 
the density profile of dark matter lead to negative exponents n, so that the density is indeed 
positive defined. 



4-. 2. Mimicking dark matter density profiles 

Many dark matter models resort to the assumption that, at the region where the galaxy rotation 
curve has flattened, both visible and dark matter might be modelled by a power law density 
profile, characterized by the outer slopes m and m', respectively (see Appendix A of Ref. [TO] 
for a thorough discussion), 
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Thus, one may establish the algebraic relation p^m ^ p^-, defining the slope ratio p = m! jm. 
The constant p^ is specific to each visible matter distribution /), and pd will be related to this 
quantity and to the previously defined pQ. Since a similar power law was obtained in the previous 
paragraph, this relation allows for a direct translation between the parameters i?0) of the non- 
minimal coupling model and the mimicked dark matter density profile quantities a', m! and p^^, 
assuming a given visible matter density profile parameters, 

1 m' ml , , 
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^^ = TT^p[p-') ^0 pS- (31) 

One may assume that the length scale a' is inherited from the visible matter density, so that 
a' ~ a. 

Clearly, the long range dominance of the mimicked dark matter over visible matter indicates 
that the density profile of the former should be shallower than that of the former, so that 
m' < m. Hence, one posits that p = m' /m < 1, thus a negative exponent n is expected. This 
range is of physical significance: since dark matter is relevant at large distances, but the scalar 
curvature increases towards the galactic core, the non- minimal coupling f2{R) should scale with 
an inverse power of the curvature. 

For concreteness, one ascertains what is the relevant range for the exponent n; to do so, one 
has to assume a particular visible matter density profile. In this study, the Hernquist distribution 
is considered [17j . 
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so that m = 4, and Eq. (|30|) leads to n = 1 — 4/m'. Bearing in mind that one aims replicating 
the dark matter distribution, it is useful to resort to two widely used candidates, namely the 
NFW profile 



with m' = 3— 7-71 = — 1/3, and the isothermal sphere profile 

with m' = 2 — )• n = —1. 

Both dark matter models are characterized by negative exponents n, as indicated above. 
The obtained values n = —1/3 and n = —1 should be taken into consideration in the following 
developments, although the general exponent n is retained. 



4-3. Background matching 

The elegant "static" solution obtained before can be fruitfully used to determine several 
characteristics of the model, as has been extensively discussed in Ref. |10) : most calculations are 
hence omitted, with the focus shifted to a comprehensive discussion of the physical significance 
of these results. 

Firstly, one notices that a negative power-law coupling /2 = (R/Rq)"' (with n < 0) does not 
admit an asymptotically fiat Minkowski space-time, since it would lead to a singular R = 
scenario. This problem is naturally avoided by noticing that the assumed coupling is valid at 
the galactic scale, and that the asymptotic behaviour should be identified with the large scale, 
cosmological background. Hence, the notion of a matching distance Too natural appears, defined 
so that the density of the mimicked dark matter approaches the order of the critical density, 
Pdm{r oo) ~ Pc (the latter being is defined as poo = Si/g / (SvrG) , with Hq ~ 70 {km/s)/Mpc the 
Hubble constant). Using Eq. (f26]l yields 
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having defined the characteristic lengthscale tq = l/^/RQ. The obtained matching distance 
which may be regarded as an indicator of the radius of the mimicked "dark" matter halo; 
for consistency, it should vanish when the non-minimal coupling Eq. (jl9|) is negligible, i.e. 
i?o — ^ if the exponent n is negative or Rq ^ oo for n > 0; this is trivially satisfied, since 
^^oo Of T^^'^ = Rq^ — ^ (excluding the constant coupling n = case). 
One may also define the characteristic density 

9 p ( c V Poo (2^GpcV 

Notice that the equation of state previously found Eq. ()28p enables a negative pressure for 
any negative exponent n; furthermore, in a cosmological setting, the non-trivial coupling 
f2{R) = (R/Rq)^ is given by (poo/po)"- Hence, the obtained relation Eq. ([36]) does hint a 
putative unified description of dark matter and dark energy, responsible for the acceleration of 
the expansion rate of the Universe. Other models of unification of dark energy and dark matter 
include the Chaplygin gas model [18] [T9[ \20[ [2T] and some field theory constructions 



4.4- Crossover between visible and dark matter dominance 

Proceeding with the analysis of the obtained "static" solution, one remarks that the model 
parameter Rq introduces both the length scale ro as well as a density scale po] the former 
should have a clear physical meaning. Indeed, rg is associated with the crossover distance Tc at 



which the density of the mimicked dark matter component becomes similar to that of visible 
matter, thus signalling the onset of the flattening of the galaxy rotation curves. By identifying 
Pdm{i'c) = p{t~c) and using the obtained profiles [271 and Eq. ([32]) . one obtains 
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with = 2GM/(? being the Schwarzschild radius of the galaxy. Hence, the crossover radius Tc 
is proportional to a sort of "geometric mean" of all the length scales present. 

4^.5. "Gradient" solution 

One ends this section by discussing the relevance of the gradient term on the r.h.s of Eq. ()23p ; 
indeed, although the "static" solution makes it vanish identically, it was obtained by neglecting 
the contribution from visible matter, since one is assuming that dark matter is dominant at 
large distances. Conversely, at smaller distances one expects the solution to Eq. ([23|) to deviate 
significantly from the "static" solution Eq. ()26p . interpolating between the latter and the classical 
identification R = SnGp. Thus, in the regime before the onset of dark matter dominance r < Tc, 
the contribution from the mentioned gradient terms may be highly relevant. 

This may be justified by inserting the GR identification R = SirGp into Eq. ()23p . and 
reading the gradient terms. To do so, one assumes that the D'Alembertian operator is given 
by its expression in a flat spacetime (in spherical coordinates), which amounts to neglecting 
relativistic corrections. This requirement translates into the condition < tq, which is easily 
satisfied by the mechanism under scrutiny, as shall be seen. 

A somewhat evolved analysis depicted in Ref. [10] indeed shows that the gradient term may 
become dominant below a distance < a < Tc, 



Notice that, at even shorter distances, relativistic corrections of the form r^/r should be 
considered. It can be shown that this is not related to the particular form of the visible matter 
density distribution (in this case, the Hernquist profile), and that assuming another form for 
p(r) would still lead to an eventual dominance of the gradient term of Eq. ()23p below a threshold 
rk- 

This concern may be attenuated by noticing that one is assuming the validity of the Newtonian 
regime: if the lengthscales a, and ro yield a value for smaller than the Schwarzschild radius, 
then the need to include relativistic corrections to the computation of the gradient term makes 
the whole approach inconsistent. As a numerical session shall show, the condition > Vg is not 
always satisfied, so that the above expression for should be considered with due care. 

Skipping over the technicalities, it may be shown that it is impossible to simultaneously 
ensure that the gradient term on the r.h.s of Eq. ()23p is negligible for all distances: it dominates 
at either r < a or r > a. This is not unexpected, given that the lengthscale a signals the shift of 
behaviour of the Hernquist profile Eq. ()32p . Hence, one is forced to conclude that the "static" 
solution Eq. ()26p is not valid after all! In particular, the discussion of the dominance of the 
gradient term in the outer region a < r < Tc is fiawed, as the crossover distance rc assumed for 
its validity. 

If the gradient term indeed dominates the r.h.s of Eq. (j23p at short ranges, it turns out that 
the corresponding "gradient" solution to this equation is the relevant quantity to assess, and 
not the "static" one already described. By the same token, it is legitimate to assume that this 
"gradient" solution will play the role of dark matter, so that it should dominate the contribution 
of visible matter at long ranges. Thus, Eq. (p3|) may be approximated by 
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One now introduces the redefined field 
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and coordinate y = rQ/r; the dominance of this "gradient" solution over the visible matter 
contribution occurs in the range r > r^, that is, y < Uk- 

In order to better understand the evolution of this "gradient" solution, one may compute the 
order of magnitude of the quantities involved at y = y^ (by recalling that, for smaller physical 
distances y > yk, the classical solution R = SnGp is valid and inserting the Hernquist profile for 
p{r)). One finds that 
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Inserting typical values for galaxies a = (1 — 10) kpc and = (10^^ — 10^^) pc (corresponding 
to galaxies of approximately (10^" — 10^^) solar masses) and recalling the validity condition for 
the Newtonian regime, rg > r^, one concludes that the r.h.s. of the above equation is less than 
unity; it turns out that the values involved are extremely small (of the order 10~^^ or lower), so 
that one may safely consider that Eq. ()39p is essentially given by 
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at large distances, i.e. small values of y — the full form of Eq. (j39p with the above redefinitions 
would be that of a generalized Emden-Fowler differential equation. Hence, the "gradient" 
solution is very well approximated by = const.; recalling Eq. (jlO]) leads to 
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that is, the "gradient" solution tracks the "static" solution! Hence, the previous results are 
partially vindicated, and the identification of the induced effect due to the non-minimal coupling 
f2{R) with a perfect fluid with equation of state Eq. ([28]) is still valid. 



4-6. TuUy-Fisher law 

The Tully-Fisher law is an empirical relationship between the intrinsic luminosity L and a power 
of the rotation velocity Voo of a spiral galaxy. An analogous relationship exists for elliptical 
galaxies, the Faber-Jackson relation. Since the luminosity is proportional to the visible mass, 
both may be written as M oc v^, with 3 < fi < 4. 

If, for simplicity, one assumes that Mdmif) oc r at large distances, thus producing a flat 
rotation curve, the familiar v'^ oc Mdm relation is recovered; resorting to the scaling Eq. ([27|l . 
this yields 

Moct;2(i-"), (44) 

so that the exponent of the Tully-Fisher or Faber-Jackson relations is related to that of the 
non- minimal coupling f2{R)- For the scenarios already discussed, one obtains M oc v^^^, for 
n = —1/3, mimicking the NFW dark matter density profile, and M oc f"^, for n = — 1, mimicking 
a dark matter isothermal sphere. 



^.7. Conservation of energy-momentum tensor 

Section [3] depicted a first attempt to describe the flattening of the galactic rotation curves as 
due to the non-conservation of the energy- momentum tensor of matter, i.e., through the non- 
geodesical motion experienced by fluid particles. Conversely, one now explores the possibility of 
mimicking "dark matter" as a consequence of the dynamics of Eq. ()23p . As explicitly stated in 
the beginning of section U one is assuming that matter follows geodesies, so that the solution 
of Eq. gives rise to the observed rotation curves. 

This assertion may now be verifled, by rewriting Eq. @, with p = for pressureless matter, 
the Lagrangian density C = —p and the adopted non-minimal coupling Xf2{R) = (R/Rq)^: 
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assuming spherical symmetry and the Newtonian approximation hrr = grr — {vr/c)'^ ~ 1- 

As discussed in Ref. [10] , a clear distinction between the "static" and the "gradient" solutions 
is that, although both are essentially the same (up to a multiplicative factor), each arises from 
widely different assumptions: the former assumes that the gradient terms on the r.h.s. of Eq. 
(I23p may be disregarded, so that the dominance of dark matter translates into the condition 
(R/Rq)^ > 1/(1 — 2n); the latter, however, is based upon the dominance of these gradient terms, 
which implies the converse relation, (R/Rq)^ < 1/(1 — 2n) < 1, for negative exponent n. 

Thus, one concludes that the issue of the covariant conservation of the energy- momentum 
tensor of matter branches between both solutions: from Eq. (|46p . one sees that the "static" 
solution would lead to 



r 



-n 



-n- 



IT 
R 



nm 



(47) 



since R oc pdm oc in the outer region of galaxies; given that n is negative, this force correctly 
points towards the galactic center, thus enabling for a flattening of the rotation curve with a 
constant asymptotic velocity. However, for small \n\ and m', this quantity would be given by 
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1, in blatant violation of the observed value Voo ~ 10^ km/ s ~ 10 ^c! Not 



only is this clearly unphysical (and would undermine the Newtonian approach taken throughout 
the text), but it also negates the whole purpose of this section, namely to assume geodesic 
motion. 

Contrariwise, the verified dominance of the "gradient" solution and related condition 
{R/RqY < 1/(1 - 2n) < 1 leads to a radial force 
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and the pathological result stemming from Eq. ()47p is avoided. 

For the already considered values of a = (1 — 10) kpc and = (10~^ — lO"'*) kpc, it suffices 
that ro > 200 kpc for the dominance of this extra force over the Newtonian one to occur only 
for distances greater than 1 Mpc. Hence, this extra force is negligible at galactic scale, and 
given the actual values for rg that will be used afterwards, much greater than this lower limit. 
Moreover, the energy-momentum tensor is conserved to a very good approximation. 



4-8. Energy conditions 

One now checks if the strong, null, weak and dominant energy conditions (SEC, NEC, WEC 
and DEC, respectively) are verified; following Ref. [23], a little algebra yields the condition 
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Clearly, both the NEC and SEC are always satisfied, for a negative power law n < 0. The 
WEC and DEC are also fulfilled for the ranges — 1 < n < and —1/2 < n < 0, respectively. 
Since the condition R > Rq is always valid at galactic scales (due to the dominance of the 
gradient term of Eq. (j23p ) for the considered models n = — 1 (isothermal sphere profile) and 
n = -1/3 (NEW profile), one has that (1 + (i?/i?o)~")~^ ~ 0; hence, the WEC and DEC are 
also satisfied for these scenarios. 

A final check is in order, namely whether the considered mechanism for the mimicking of dark 
matter does not give rise to the so-called Dolgov-Kawasaki instabilities, whereas small curvature 
perturbations could arise and expand uncontrollably [21] ■ Refs. |23t I25j allow one to write the 
mass scale tudk for these as 
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which are given approximately by 



Since rnj-fj^ is negative, the possibility of instabilities at a length scale r^K is worrisome, with 

(R \ -'-^"/^ I — 
it) \ ~'^o- (53) 
R J \J Po 

At large distances, where the dark matter-like component dominates and the scalar curvature is 
greater than the matter density, one finds that this typical length scale r^K grows linearly for the 
m' = 2 isothermal sphere scenario, and with a power-law r^/^ for the m' = 3/2, corresponding 
to the NEW density profile. 

It should be noted that the generation of these instabilities assumes an initial small 
perturbation in a constant background curvature; this is best suited for an analysis in cosmology, 
but should be approached with due care in the present context, where there is a spatial variation 
of R{r). Indeed, one must instead resort to the characteristic lengthscale over which one may 
assume a constant curvature, given by L = —R[r)/R'[r) = r/m'. One may compare this 



constant curvature domain L with the instabihty length scale r£,K to ascertain if the assumed 
perturbative expansion, as embodied in the inequality r£,K < L, is valid. Equating both 
quantities yields 

Q-<,-»«'^(2*.)^j^^\/^. (54, 

This relationship immediately shows that the perturbative expansion is invalid if n = — 1, when 
dark matter assumes a isothermal sphere distribution; if it takes the form of a NFW profile, 
n = —1/3, one finds that the perturbative expansion is valid only below distances r < r^, 
where relativistic corrections become relevant. Hence, the Dolgov-Kawasaki instabilities are not 
physically meaningful in the context of this work. 



5. Mimicking the dark matter density profile of spherical galaxies 

This section aims to apply the developed mechanism for the mimicking of dark matter to fit 
numerically several rotation curves of galaxies of EO type — these are chosen since they exhibit 
an approximate spherical symmetry adequate for the results obtained in this work. Seven galaxy 
rotation curves, are used, with the division into visible and dark matter components as reported 
in Refs. [Ml EZl EEl [29] : NGC 2434, NGC 5846, NGC 6703, NGC 7145, NGC 7192, NGC 7507 
and NGC 7626 galaxies. 

As stated in the preceding section, particular interest is assigned for n = —1 and n = —1/3, 
which correspond to the isothermal sphere or NFW dark matter density profiles. Furthermore, 
one may recall that the power-law profile f2{R) = (R/Rq)^ was discussed in the context of a 
possible series expansion of a more general non-minimal coupling; hence, there is no reason to 
select a priori either case n = — lorn = — 1/3: instead, one can successfully take advantage of 
a composite coupling given by 



(55) 



so that both profiles may be added in order to describe the dark matter density profile and the 
related galaxy rotation curve. 

The proposed non-minimal coupling leads to a more complex expression for the trace of the 
field equations 
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One can no longer resort to the identification of a "static" solution, as before, since the non- 
linearity of the above equation indicates that the solution for this composite scenario is not 
given by the mere addition of individual solutions; by the same token, no clear equation of state 
or identification of the mimicked dark matter profile with a perfect fluid arises. Nevertheless, 
one expects that the main feature of the model are kept, namely the overall mimicking of dark 
matter density proflles according to Eq. (j27p and the dominance of the gradient term on the 
r.h.s. of the above equation. 

This stated, one proceeds as follows: each visible matter density profile p(r) is derived from 
the visible component of the corresponding galaxy's rotation curve; two added Hernquist profiles 
are used, one for the extended gas and another for the core matter distribution. This density 
profile is then used as an input to Eq. ()56p . and a best fit scenario for the dark matter component 
of the rotation curve is obtained by varying the parameters ri = 1/y/Ri and rs = l/^/R^. 
Although these should be universal quantities, one resorts to individual fits of each galaxy. 
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Table 1. Best fit values for the cliaracteristic lengtliscales ri and for the composite and 
separate fits of the galaxy rotation curves with the n = —1 isothermal sphere and n = —1/3 
NFW mimicked "dark matter" scenarios, together with background matching distances i 
(n = —1) and Too 3 (^^ = —1/3) for the composite non-minimal coupling. The limit rj = oo for 
the composite scenario indicates that the corresponding scale Ri = . The units used are ri 
(Gpc), rs (10^ Gpc), Too 1 and r^o 3 (kpc). 



allowing for a posterior discussion of order of magnitude, deviation between obtained values, 
any abnormal case or possible trends, etc. 

The results derived from the numerical solution of Eq. (|56p are depicted in Table [H and 
shown in FigH} as can be seen, the composite non-minimal coupling Eq. (|55|) provides close 
fits for all rotation curves in the outer region — of much higher quality than those obtainable if 
one assumed individual couplings f2{R) = R/Ri or f2{R) = ^/R/Rs (see additional graphs in 
Ref. [To]). For some galaxies, the overall quality of fit is disturbed by a discrepancy in the inner 
galactic region. This could result from an uncertainty in the original derivation of the rotation 
curves, or perhaps due to the deviation from purely spherical symmetry. 

From a fundamental standpoint, this deviation could also indicate that the non-minimal 
coupling should include an additional term that acquires relevance in a higher density 
environment: recalling an earlier work in the high density context of the Sun, one could admit a 
linear addition 6f2{R) oc R, for instance. Of course, the considered model could be degenerate 
to the approximation of the pure curvature term given by the GR prescription fi{R) = 2kR, 
and a more complex picture involving both non-trivial fi{R) and f2{R) couplings is in order. 

Although details are here suppressed, for brevity, after obtaining the solution to Eq. ([56]) it 
is straightforward to compute the gradient term on its r.h.s.: as discussed before, this is indeed 
dominant, so that the mimicked dark matter density profile is almost completely overlapped 
with it. 

5.1. Characteristic density, background matching and mimicked "dark matter" inter- dominance 
The selected galaxies have their rotation curves fitted by characteristic lengthscales of the order 
ri ~ 10 Gpc and r^ ~ (10^ — 10^) Gpc. Although the analytic results derived in the preceding 
sections do not hold strictly in the assumed composite coupling scenario of Eq. ([55]) . Eq. (j35]) 
should nevertheless provide a quantitative figure for the background matching distances. 
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Figure 1. Left: Observed rotation curve (dashed full), decomposed into visible (dotted) and 
dark matter (dashed grey) contributions [26], superimposed with the mimicked dark matter 
profile (full grey) arising from the composite non-minimal coupling and resulting full rotation 
curve (full). Right: Log-Log profile of the visible matter density (dotted), mimicked "dark 
matter" contribution (dashed) and the sum of both components (full). 
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n = — 1 



after inserting the Hubble constant = 70.8 km s M-pc . Notice that this expression 
was obtained with the assumption of a single Hernquist profile for the visible matter density, 
and not for the two component model used in the current fitting session; using the dominant 
contribution stemming from the extended gas distribution, one obtains the values displayed in 
the last columns of Table [TJ 

Since may be regarded as signaling the radius of the mimicked dark matter haloes. Table [T] 
shows that the n = — 1/3 NFW haloes are always smaller than those with the n = — 1 isothermal 
sphere profile — with the exception of NGC 2434 that, however, stands out from the set since it 
is best fitted by a single NFW "dark matter" component. Also, r^Q i is always larger than the 
range of the rotation curves, so that the isothermal sphere "dark matter" extends further than 
the observationally inferred; on the contrary, the value of 3 goes from about 0.57-L (NGC 
7145) to 1.73L (NGC 7507), L being the full distance considered for the rotation curves. 

This indicates that the isothermal sphere halo is more significant that its NFW equivalent: 
even the exceptional case of NGC 2434, where no isothermal sphere component appears in the 
best fit scenario, yields a matching distance roo 3 = 33.1 kpc — only about twice the endpoint 
of the corresponding rotation curve, L ~ 17 kpc. 

This is confirmed by a simple algebraic computation: since the NFW profile falls as r~^, while 
the isothermal sphere one displays a shallower long-range dependence, and the respective 
characteristic densities obey p3 ^ pi , it is easy to conclude that the isothermal sphere eventually 
becomes the dominant dark matter contribution. Indeed, equating both contributions yields 



so that this dominance occurs at distances larger than r = a/10 ~ 1 kpc. Thus, the considered 
galactic rotation curves are found to be "mostly flat" . 

5.2. Cosmological relevance 

As discussed already, the results obtained above are naturally tuned towards the corresponding 
astrophysical setting, if one assumes that the chosen form for the (composite) non-minimal 
coupling represents the leading terms of a series expansion of a more evolved expression for 
/2(-R), suitable in the mentioned range for the scalar curvature R. Hence, there is nothing 
mandatory concerning a possible cosmological relevance of the composite coupling Eq. (j55p . with 
the obtained values for Ri and R3: it might be the case that a suitable mechanism applicable 
in cosmology corresponds to another term of the hypothetical series expansion of f2{R)- 

Nevertheless, it is natural to assess if the considered cases n = — 1 and n = — 1/3 give rise to 
implications for cosmology, since this would present a natural candidate for a unification model 
of dark matter and dark energy. From the discussion of paragraph 14.31 one may write the non- 
minimal coupling Eq. ()55p . in a cosmological context, as 



having used R = 'i/r'jj, with rn = c/H = 4.2 Gpc the Hubble radius, and inserting the best fit 
orders of magnitude ri ~ 10 Gpc and ~ (10^ — 10^) Gpc. 
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The small value of the n = —1/3 non-mmuTial coupling indicates that this should not be 
significant at a cosmological level; the higher value of the n = — 1 power-law coupling (when 
compared to unity) hints that a simple Ri/R coupling could prove interesting for cosmology. 

5.3. Universality 

The obtained values for ri and r^, shown in Table [H display an undesired variation between 
galaxies, since these parameters should be universal. Indeed, ri averages fi = 21.5 Gpc with a 
standard deviation ai = 10.0 Gpc, while rs presents = 1.69 x 10® Gpc and o"3 = 1.72 x 10® Gpc. 
This deviation from universality is not unexpected, and several causes could be hinted as culprits: 
although small, the existing assymetry in the selected type EO galaxies due to the 1 — 6/a < 0.1 
semi-axis deviation from pure sphericity could translate into the difference in obtained values for 
ri and r^. This shift from the assumed symmetry could be enhanced by localized features and 
inhomogeneities, which could also act as "seeds" for the dynamical generation of the mimicked 
"dark matter" component. Furthermore, the reported dominance of the gradient term in the 
r.h.s. of Eq. ()56p could boost this effect, if these perturbations vary significantly on a short 
lengthscale. 

By the same token, the relevance of the gradient term could lead to the deviation from 
universality if the chosen profiles (Hernquist for visible matter, NEW and isothermal sphere 
for mimicked dark matter) prove not to be so suitable as hoped: this mismatch between the 
real density profiles and the fitted ones does not need to be large, if a large deviation of its 
first and second derivatives lies in very localized regions: this could act as a sort of "fitting 
inhomogeneities" with the same effect as of those discussed above. 

Aside from the symmetry and fitting functions selected, more fundamental issues could be 
behind the variation of the obtained values for ri and r-^: in fact, the non-minimal coupling Eq. 
(|55p could be too simplistic, and a more complex model could provide a more universal fit of 
the galaxy rotation curves. 

Also (as discussed with respect to the relatively less successful fits of the inner galactic 
regions), the simplifying assumption of a linear curvature term fi{R) = 2kR might also 
contribute to this deviation from universality, since one knows that power-law terms fi{R) oc R"" 
can also induce a dynamical "dark matter" |8j. 



Finally, a more interesting idea could lie behind the obtained discrepancy, namely that the 
visible matter Lagrangian density itself may be different from the assumed perfect fiuid form, 
= —p [13(, |14J: this functional dependence is inadequate in the context of the model here 
addressed, and a more adequate form Cm = Cm{p) should be considered — perhaps even 
including other variables related to the thermodynamical description of the matter distribution 
(see Ref. [30] for an in-depth discussion based upon a similar consideration). 

6. Conclusions 

In this work the possibility of obtaining a solution to the dark matter puzzle, embodied by the 
fiattening of galaxy rotation curves, was approached by resorting to the main phenomenological 
implications of models possessing a non-minimal coupling of matter to curvature, following 
results reported in Ref. [TO]. 

As a first attempt, one first examined the non-conservation of the energy- momentum tensor 
of matter and the implied deviation from geodesic motion, and ascertained what coupling f2{R) 
should be so that the derived extra force would lead to the reported fiattening of the rotation 
curves. This requires a logarithmic coupling of the form Xf2{R) = —v^/m\og[R/R^,) (where 
m is the outer slope of the visible matter density p), which can be approximated by a simpler 
power-law Xf2{R) ~ (i?*/i?)°, with the asymptotic velocity given by = am. However, 
this solution yields an almost universal Voo, prompting for another possible mechanism for the 
dynamical mimicking of "dark matter" . 



This need led to the second, more evolved approach, which is rooted empirically in the 
phenomenological Tully-Fisher relation: one may instead assume that geodesical motion is 
preserved, V^T^,y = (a condition proved to be self-consistent), but that the metric itself 
is perturbed: the mimicked dark matter density is then given by the difference R/2k, — p. By 
resorting to a power-law coupling with matter f2{R) = [R/Rq)^ (with a negative exponent n 
yielding the desired effect at low curvatures and long range). 

The proposed power-law directly leads to a "dark matter" density that depends on a power 
of the visible matter density, thus accounting for the Tully-Fisher law in a natural way. The 
obtained "dark matter" component has a negative pressure, as commonly found in cosmology 
in dark energy models: this hints the possibility that the non-minimal coupling model might 
unify dark matter and dark energy. 

Given their relevance in the literature, two different scenarios were considered: the NFW and 
isothermal sphere dark matter profiles (n = — 1/3 and n = — 1, respectively). Although separate 
fits of the selected galaxy rotation curves to each coupling do not yield satisfactory results, a 
composite coupling of both power-laws produced a much improved adjustment. 

The characteristic lengthscales ri and rs were taken as fitting parameters for each individual 
galaxy, yielding the order of magnitudes ri ~ 10 Gpc and ~ 10^ Gpc. This enabled the 
computation of the cosmological background matching distances for each galaxy (which depend 
on their characteristic lengthscale a) and the obtained astrophysical range showed that the 
n = — 1 isothermal sphere "dark matter" halo dominates the n = —1/3 NFW component. 
Furthermore, the n = — 1 scenario was shown to have possible relevance in a cosmological 
context, as ri ~ r//, the latter being the Hubble radius — again hinting at a possible unification 
of the dark components of the Universe. 

The lack of the desirable universality in the model parameters Ri and it's was duely noticed, 
although these quantities are characterized by the same order of magnitude for the best fits 
obtained; several possible causes for this variation were put forward, from deviation from 
spherical symmetry to a more complex form for the non-minimal coupling between curvature 
and matter f2{R) or the curvature term fi{R), or the need for a more evolved Lagrangian 
description of the latter. 

As a final remark, one concludes that the rich phenomenology that springs from the model 
with a non-minimal coupling between matter and curvature enables a direct and elegant 
alternative to standard dark matter scenarios and many modifications to GR, which usually 
resort to extensive use of additional fields and other ad-hoc features. As an example of the 
latter, one points out the MOdified Newtonian Dynamics (MOND) hypothesis, which is by itself 
purely phenomenological, and whose underlying Tensor- Vector-Scalar theory is based upon an 
extensive paraphernalia of vector and scalar fields [311 [321 ESI El] (see Ref. [|35j for a critical 
assessment). As a final remark, we point out that the considered model may also be translated 
into a multi-scalar theory [12] . with two scalar fields given by 
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with dynamics driven by a potential 
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